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Introduction 

This is the final technical report for NASA grant NAG 8—431 "Limits on 
Coronal Material in Normal Galaxies." 

The original objective of this project was to determine whether useful 
information about the diffuse thermal emission from other galaxies could be 
obtained from Einstein IPC images. This goal was reached after considerable 
effort to understand the IPC background and how to correct for it. The work 
has resulted in one talk and two published papers. 

The results obtained from an analysis of the Ml 01 images could be used 
to rule out for Ml 01 some models of diffuse thermal emission that had been 
proposed for our galaxy. More extensive analysis in terms of models of 
supernova remnant evolution in a galactic disk shows that additional 
constraints can be placed on the effective density of the galactic disk in 
regions where supernovie occur. This work was done in collaboration with 
Don Cox, and has been published in the second papar (a copy of this paper in 
its final form is attached to this report). 

The most serious limitation on the conclusions reached was the size of 
the statistical sample: one galaxy. There were also several other 
interesting results that were statistically marginal or just out of reach. 
Unfortunately, the Ml 01 images are by far the best data of this type taken 
by Einstein , and little more can be done with this data base. ROSAT and 
AXAF will have far better capabilities for these purposes, however, and in 
the course of this Einstein data analysis project we have learned both the 
range of interesting things that can be done with this kind of observations 
and the most efficient way to use the new satellites to do them. 
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ABSTRACT 

Measurements of the X-ray surface brightness of a face on 
disk galaxy M101, have previously been used to place upper limits 
on the power radiated by a hot corona. Such analysis constrains 
the effective density of the disk; either it must be so low that 
the remnants drive a fast hot wind (low radiated power) or so high 
that the remnant temperature at overlap is low (low X-ray 
power). These X-ray measurements are here used to constrain the 
properties of the population of supernova remnants evolving in the 
disk. This adds a further constraint since young remnants 
evolving in higher density radiate more of their energy in X-rays, 
whether or not they eventually overlap to generate a hot corona. 

The strength of this second limit depends strongly on the 
density history of the remnants (e.g., evaporative versus 
nonevaporative evolution) and on the assumed supernova rate. For 
evaporative evolution the analysis rules out the McKee and 
Ostriker ISM model in particular and evaporative evolution in 
general unless the supernova rate is at least several times lower 
than current expectations. For standard Sedov evolutions, the 
density limit marginally admits evolutions in 0.2 cm -3 , a popular 
alternative to the McKee and Ostriker model. 
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I. INTRODUCTION 

Two attempts to use the Einstein satellite IPC to measure the 
X-ray emission of fountains, coronae, or hot interstellar gas of 
noncluster disk galaxies have made no clear detection. Bregman 
and Glassgold (1982) studied two edge on spirals (NGC 3628 and NGC 
4244), a very sensitive test for coronae and fountains, while 
McCammon and Sanders (1984) tried the face on spiral M101, for 
which X-ray emitting material dispersed through the interstellar 
medium could potentially be seen. In both cases it was found that 
the upper limits on the X-ray emission were some two orders of 
magnitude lower than the supernova power estimated to be 
available. Essentially the same conclusion applies to a local 
region of our own galaxy where the measured surface brightness of 
the soft X-ray background is a similarly small fraction of the 
mean surface density of supernova power (e.g. Cox, 1981, 1983). 

These results restrict the manner in which supernova energy 
is dissipated, thus limiting the range of tenable models for 
supernova remnant evolution in the interstellar medium. At a 
given supernova rate, the model parameter most directly 
constrained is the degree to which remnants, on average, dilute 
their energy with interstellar material before overlapping one 
another or entering the radiative phase. 

The options available prior to these restrictions seem to 


have been: 
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1) Little dilution (less than ~100 M # per supernova, very low 
density ISM and little evaporation); SN energy drives a hot low 
density wind from the galaxy. Owing to the low density. X-ray 
emission is weak. 

2) Intermediate dilution (a few hundred solar masses per 
supernova); SN generate an active fountain or halo which radiates 
the entire SN power in X-rays and EUV. The fraction of the power 
radiated in X-rays decreases with increasing dilution (more mass 
dilution implies lower coronal temperature). 

3) Critical dilution (roughly 3000 M # per SN, acquired late 
at low effective density); individual remnants radiate the bulk of 
their energies shortly before overlap with one another, providing 
a large filling fraction of hot gas in the disk. The coronal 
component temperature is so low that it emits few X-rays. This is 
the general form of the McKee and Ostriker model (1977, hereafter 
MO). Or finally, 

4) Large dilution (high density ISM or much evaporation) 
individual remnants become radiative while still relatively small, 
forming large expanding shells that snowplow through the medium 
until encountering one another long after the transition to the 
radiative phase. 

One might suppose that for dilution even greater than the 
critical value. X-ray detectors would be insensitive probes. For 
observations of a face on spiral galaxy, however, such is not the 
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case. A detector can look into the disk and record the integrated 
surface brightness of the population of individual remnants. Even 
at critical dilution, X— ray emission by young hotter remnants can 
be significant. For greater ISM density or higher thermal 
evaporation from clouds (in either case greater density of hot 
material within remnants, greater dilution), the X-ray emission of 
the individual remnants in the population is higher. 

For a given supernova rate, the X-ray surface brightness can 
thus be used to set limits on both coronal and remnant population 
properties, the former constraining the intermediate dilution 
possibilities (setting a lower limit on dilution in that regime), 
the latter constraining higher dilutions, setting an upper 
limit. For an important subset of possibilities (the assumption 
of evaporative remnants), these two limits on M101 turn out to be 
mutually exclusive, implying either that the supernova rate is 
lower than expected or hot remnants evolve in a less X-ray 
emissive fashion. 

The average galactic distribution of emission measure (versus 
temperature) from its population of individual remnants is 
calculated in Section II. Section III reviews the analysis of the 
Einstein IPC measurements of M101 by McCammon and Sanders (1984) 
and presents the upper limit on the remnant population properties. 
The latter is derived by restricting the calculated X-ray 
brightness, from the emission measure distribution of Section II, 
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to be less than the observational limit. The surface brightness 
depends on four parameters (supernova power, characteristic 
remnant density in the X-ray emitting regime, density evolution of 
remnants with age, and overburden of X-ray absorbing material in 
M101). For reasonable estimates of the supernova power and 
absorption, the limit is expressed as a maximum remnant internal 
density in the X-ray emitting regime (specifically at 10 6 K) , as a 
function of the assumed mode of remnant evolution. 

For remnant evolution with constant or increasing average 
density, the extreme upper limits on remnant internal densities 
are comparable to the average interstellar density in the Milky 
Way, certainly exceeding the expected density within large 
supernova remnants. For some reasonable parameter choices, 
however, the constant average dehsity case marginally threatens to 
exclude an important possibility. For evaporative remnants, whose 
mean internal density decreases with radius, the density limit is 
more stringent. The broader implications of the latter are 
considered in Section IV which studies both the remnant population 
and coronal emission limits, deriving the properties of an MO 
state which is just barely consistent with the X-ray limit. The 
results are summarized in Figure 2. Section V presents an 
overview of the results and discusses the importance of future 
observations on studies of this type. 
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II. THE EMISSION MEASURE DISTRIBUTION OF A POPULATION OF REMNANTS 

Although an isolated supernova remnant may radiate the bulk 
of its energy only after evolving to a rather low temperature (< 4 
x 10**K) , it nevertheless passes through a phase at X-ray emitting 
temperatures and radiates some of its energy there. As a result, 
a steady state population of such remnants will always have some 
members contributing to a high temperature tail of the emission 
measure distribution function. 

For this study we will assume that individual remnants mature 
in an isotropic environment and radiate as though they were 
completely in collisional equilibrium at the post shock 
temperature. The effect of this assumption can be assessed from 
the results of Cox and Anderson (1982) and Hamilton, Sarazin, and 
Chevalier (1983). Using the equilibrium rates slightly 
underestimates the X-ray emissivity. Assuming that the emission 
spectra can be represented by the post shock temperature will 
ignore emission from the hot interiors of remnants which have 
slowed too much to generate X-rays at the edge. More detailed 
calculations would therefore provide more restrictive upper limits 
on the high temperature tail of the emission measure distribution 
and the effective density of X-ray emitting remnants. 

The hot mass within a supernova remnant during the adiabatic 
phase will be proportional to radius cubed only in the simplest 


case. Generalizing to 
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M « R 3+b (1) 

s 

where b ■ -5/3 is the MO thermal conduction value, b - 0 for a 
homogeneous medium, and b is positive for evolution in a cavity 
where density increases gradually with radius, we then have 

T « v 2 « E /M, R « t 2/(5+b) , and v ** [2/(5+b)] R /t. 

S S 0 8 S S 

It is straightforward to express all of the evolving parameters of 
a remnant as power laws of the post shock temperature. An 
important example is the distribution of remnant ages over 
temperature : 


dt = 5+b t 

dT = 2(3+b) T. 

The luminosity of a single remnant can be written (e.g., Cox 
1985a) 


£ m 8(b, q) N n A L(T g ) 

or (3) 

£ - B’(b,q) n 2 V L(T ) 

s 

where V is the remnant volume, N = 3Vn A /(3+b), n A is the effective 
preshock ambient density (the actual value except in evaporative 



9 


models), n * 3n^/(3+b) is the average density of hot gas within 

the remnant, the parameter q arises from the assumption that the 

cooling function L(T) is proportional to T -< ^ and 8 and 8' = (3+b)8/3 

are effective compaction parameters (see Cox and Franco 1981, or 

Cox 1985a for examples). Dividing £ by L(T g ) provides the 

2 

effective value of / n dV at T . 

1 e s 

In a disk galaxy of area Ag and supernova rate (t sn ) *, the 
number of remnants in the post shock temperature range T to T+dT 
per unit area of the disk is 


dv 


A g t sn 


dt 

dT 


dT 


(4) 


so that the distribution function over temperature of surface 
averaged emission measure is 


d (EM) m £(T) 1 5+b t_ 

dT = L(T) A g t SN 2 ( 3+b ) T 

2 

= 3(5+b) 8 n A (T) V(T)t(T) 

2(3+b) 2 A g t SN T 

Choosing a particular temperature T Q at which to normalize 
these results, they can be written entirely as functions of 
conditions when T g *= T Q (all shown with subscript zero). 


multiplied by a power law in T g /T 0 : 
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d(EM) 

d(T/T ) 
o 


_J r 3(5+b)g v n2 

A g T SN 2(3+b) 2 ° Ao 


17+7b 

tj (^) 2<3tb) 


( 6 ) 


By using the numerical integrals in Cox and Franco 1981 (see also 
Cox (1985a) the square bracket for b > 0 can be approximated 


3 < 5 + b >g v n 2 t 
2(3+b) 2 o Ao o 


1.2 x 1 0 72 cm -3 s G(b) E 51 4 / 3 ff Q 2 / 3 


• (10 6 K/T o ) H / 6 


(7) 


where 


G(b) 


JL9 

36 


(1 + 0.239b) 


4/3 


( 8 ) 


defined such that G(0) * 1 . The function 8(b,q) was evaluated 
for q ** -1/2 by Cox (1985a) from the Cox and Franco (1981) 
integrals to be (for positive b) : 



1 2 ( 3+b ) 
l9+8b * 


(9) 


For b “ -5/3, the evaporative case, 0' * 2 (implying 8 “ 4.5) was 
advocated by Cox (1985a), because the above formula from fitting 
positive b results suggests an absurdly small compaction parameter 
at b « -5/3. Choosing the normalization to be T Q * 10 b K to place 
it in the X-ray regime, and expressing the supernova rate per unit 
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area relative to 1 SN per 30 years within a disk of radius 15kpc, 
we finally have 


6 17+7b 

d(EM) „ ^ , n ~2 -6 x „„ ^ rlO 2(3+b) 

— £ — - (6,3 x 10 cm pc) G(b) [ — — J 

d(T/10 K) 

. (ISkECj 2 E5 i 4/ 3 - 2/3 . (10) 


For fixed n Q , this distribution weights high temperatures more 
heavily for negative values of b because the remnants are then 
denser when hot. Notice that n^ Q is specifically the effective 
preshock density when the post shock temperature is 10^K, 
while n Q is the average density within the remnants at that 
temperature. Depending on b, the densities can be either higher 
or lower at earlier epochs. A population of remnants with a 
distribution function for n Q simply introduces <n 0 2 / 3 >, while a 
distribution of b values would be more complicated. 

An important caveat regarding the above approximation is that 
remnants have been assumed to be evolving in an isotropic 
environment so that their luminosities and radial evolution are 
coupled via n^(R g ). The inclusion of significant inhomogenity in 
the density distribution on a scale comparable to the remnant 
size, like that observed, for example, in the Cygnus Loop, 
complicates the analysis considerably, making the density values 
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inferred in this study measures of something intermediate between 
the most pervasive density and the density of the brightest X-ray 
emitting features. Such inhomogeneities can be regarded as 
enhancing the effective value of the compaction parameter g'. 

Since this parameter provides the conversion between rms and 
average density, using the value appropriate to isotropic remnants 
is conservative in evaluating the upper limit to the average 
density within remnants. 
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III. COMPARISON WITH Ml 01 

The Einstein satellite IPC measurements of M101 were analyzed 
by McCammon and Sanders (1984, McS) In concentric rings of outer 
radius 5i arcmin, where i-1,5. The visible galaxy is prominent in 
rings 1 and 2. The three historical supernovae occurred In rings 
2 and 3. The Holmberg radius lies approximately at the outer 
boundary of ring 3. Data in the vicinity of identifiable sources 
were omitted. 

The results of these observations were that: 

(1) Rings 2 through 5 have essentially identical soft X-ray 
surface brightness, all slightly lower than anticipated from the 
rocket measurements of the 8° average soft X-ray background in 
that direction. The background is unusually bright in this 
region, but there is no indication that it requires other than the 
usual origin in the solar vicinity. The surface average emission 
measure of M101 hot gas in these outer rings is certainly less 
than that of the Local Bubble. 

(2) Ring 1 has a detectable excess count rate but, comparing 
IPC exposures separated by 6 months, some of the soft X-ray flux 
is clearly varying in time. It is not known how much of the 
remaining ring 1 flux should be attributed to non- SNR sources. 

The observed count rates and their lo uncertainties are in 
Table 1 . By considering rings 4 and 5 as background (non X-ray 
plus foreground emission) measurements, upper limits can be set on 
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the SNR induced count rate for rings 1 through 3. Roughly 
speaking (depending on channel), the 3a limits are then 80%, 10%, 
and 10% respectively, of the total observed rates in rings 2 
through 5. 

The emission measure distribution in equation (10), combined 
with Raymond and Smith (1979) equilibrium spectra and absorption 
due to assumed intervening material, is folded through the IPC 
response function to calculate anticipated count rates in each IPC 
channel. This must be done separately for each interesting value 
of b (choice of SN evolution mode) and (assumed column density 
of absorbing material). The calculation performed was normalized 
to 

B $ (il^E£) 2 (iSHi) E 4/3 72/3. (11) 

E 2 t SN 51 

g 

Notice that B is proportional to the supernova power per unit area 
multiplied by the cube root of the Chevalier scaling parameter 
En. The integral over the temperature distribution of equation 
(10) can be extended to T ** 0 because remnant cooling, at 
densities of interest, does not set in to alter or truncate the 
distribution until shock temperatures below those to which the IPC 
is sensitive. 

By weighting the channels of greater count rate more heavily 
and imposing overall consistency at the 3a level with the set of 
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rates and their standard deviations, one finds the 3a upper limit 

to the normalization parameter B. Two cases were considered, ring 

1 alone (minus the average background of rings 4 and 5) and an 

area weighted average of rings 2 and 3 (again minus 4, 5). The 

results are shown in Figure 1 where the hydrogen column density is 

in units of 10^0 cm - 2 an< j the indicated value is in addition to 
20 -2 

the 1.1 x 10 cm contribution of the Milky Way in the M101 
direction. As expected from equation (10), the limit on B is 
tighter for smaller b. In addition, the limit on B for rings 2 
and 3 is consistently a factor of about nine more stringent than 
for ring 1, scaling with the respective X-ray limits. Owing to 
the unknown contribution from non-diffuse sources, ring 1 will not 
be considered further. 

If remnants evolve in very 'low densities so that they achieve 

large sizes, then a small Njj overburden is expected, at least for 

20 —2 

the near side. (Our Local Bubble has ~1.2 x 10 cm .) 

Conversely, remnants evolving in higher density remain small and 

20 —2 

more likely to experience half of the 6 x 10 cm full disk 
column density characteristic of Ml 01 as a whole. 

• Figure 1 should serve as a constant reminder of the tenuous 
nature of the conclusions about to be presented. Continuing in a 
conservative vein, we discuss only the upper limit provided by the 
Njj » 3 x 10^® cra - ^ results. 
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From Figure 1 we infer that for evaporative evolution, B < 
0.01; for evolution at constant average density, B < 0.08; while 
for evolution with mean density increasing as radius cubed, B < 
0.2. Writing the normalized supernova power per unit area 


P - E (i|!sE £) 2 (iS=£) 
g t sn 


( 12 ) 


such that 


P(E 


- 2 ^ 1/3 
o 


51 n - > 


(13) 


we can then write the effective remnant density at the epoch T s = 
10 6 K 

n Q * = <n o 2 / 3 >3/2..= ( B / P ) 3/ 2 ^ 1 / 2 (14) 

Uncertainty in the supernova explosion energy introduces at least 
a factor of four uncertainty in density limits. Once again taking 
the conservative position with £ 5 } >0.5, T gN < 100 years per 
15kpc radius area, we have P > 0.15, and 

n ' < 24 B 3/2 . (15a) 

o ~ 

For b = -5/3, 0, and 3 the extreme upper limits on n ' are 

__o 

then 0.03, 0.6, and 2 cm , respectively. Unless the supernova 
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rate is especially low, these are the maximum allowable average 
densities of hot gas within remnants whose shock temperature is 
10 6 K. 

From a study of this type one is interested in learning 
whether remnants might conceivably evolve in an intercloud medium 
of density 0.1 to 0.3 cm -3 or whether they necessarily evolve in a 
much lower density component , with or without appreciable thermal 
evaporation of embedded clouds. For b >0, it is clear that 
limits posed by the Einstein observations are too weak to disallow 
any interesting possibility for E 51 = 0.5 and P = 0.15. In fact, 
the b=0 standard Sedov case has the more general result 


o' < °:° 23 r 3 ( ( ->-t 3/2 

o E * ‘15kpc‘ ‘30 yrs * 


(b= 0 ) 


(15b) 


so that lnslstance on E^ = 1 and contributing supernova rate per 
unit area comparable to that in the Milky Way does endanger the 
possibility of a moderate density intercloud component. With 
improved observations discussed in Section V, a definitive 
statement may be possible. In the meantime, the limit for the 
evaporative case is already so strong that only low density 
remnants are compatible. The limits from individual remnants and 
hot coronae are combined in the next section, showing the severity 
of the problem with this mode. 
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IV. EVAPORATIVE EVOLUTION 

The X-ray surface brightness of a galactic disk whose 
supernova remnants undergo evaporative evolution depends on three 
parameters: the supernova power per unit area (P), the effective 
density of remnants at some fiducial temperature (n Q ) > and the 
column density of absorbing material. Having taken half the total 
column density as representative of the latter, we now pursue the 
regime in P, n Q space allowed by B < 10“ 2 . In addition to 
criteria separating allowed from disallowed regions, we shall 
present those distinguishing intermediate from high dilution 
regimes. The plane and its subdivisions are represented in Figure 
2 . 

Approximate remnant properties as functions of E 0 , b, T g , 
and n were presented in Cox (1985a). For b = -5/3, the needed 
results include the extrapolated adiabatic phase temperature at 
the cooling epoch: 


T„ = 0.88 x 10 6 K (E C1 n 2 
c 51c 


(16a) 


0.82 x 10 6 K (E 51 n 2 ) 2/9 , 


and the corresponding radiated fraction of the original energy 


0 . 68 , 


(16b) 
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where n c is the average remnant density at that epoch. The 
general relation between density and temperature is 


n 


% t 6 5/4 - ”c <V t 6,c> 5/4 ' 


The remnant radius as a function of n and Tg ** T g /10^K is 

E s i i / o 

R g = 22.8pc [ _ (17) 

n 6 


while the areal porosity factor is 

Q A = f irR g 2 dt/(*R g 2 t sn ) = 1.75 x 10~ 3 - P - (18) 


The latter is related to the radiated fraction by 


Qa= Q A ,c [ f /°.6S] 11/6 . (19) 

We continue the notation using B and P of equations (11) and (12). 

The parameters of a population that behaved approximately in 
accordance with the McKee and Ostriker (1977) ISM model (critical 
dilution) can be inferred by forcing Q A to be 1/2 at the cooling 
epoch. They include: the average density, temperature, and radius 
at that epoch 


- 0.022 cm 


-3 


n 


3720 


(TFT--) 


7/10 


J 51 


A,c 


( 20 ) 



20 


T 


6 ,c 


0.29 


*■51 


1/10 




1/5 


( 21 ) 


s ,c 


123pc E 51 7/2 ° ( 


2Q 


A,Cn3/10. 

J » 


( 22 ) 


the average density and radius at 10 b K 


n 

o 


0.10 cm 

Z 11 MO 

E 51 




9/20 


A,c 


(23) 


R so = 49 pc 


"51 


17/40 


r 2Q A,c^3/20. 
I P J ’ 


(24) 


and the brightness factor for the population 


B = 0.22 


F 3 / 20 P 

E 51 P 


^^A.c 5 


13/10 
3710 * 


(25) 


In Figure 2, the locus of MO states is provided by equation (23) 
with Q a c = 0.5. 


Although the densities, temperatures, and radii seem quite 
reasonable for P ~ 1, it is clear that the X-ray emission would in 
that case be far greater than the upper limit for Ml 01. It is in 
this sense that the X-ray limit is restrictive of the critical 
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dilution possibilities. Consistency is possible only if the 
supernova rate is small. 

The MO state providing a brightness factor B requires 

P E 51 3 / 26 - 0.09 (100B) 10 / 13 (26) 


and 


T SN " 330 5™ E 51 


29/26 ^ 15kpc j2 


1 


g (100B) 


wn‘ 


(27) 


Thus, since 100 B < 1, our upper limit on the X-ray emission of 
this MO state population is consistent only with a supernova rate 
per unit area four to eleven times smaller (depending on E^) than 
that commonly assumed for the Milky Way. In Figure 2, inter- 
section B is provided by equation (26) with 100 B = 1. 

The evaporative remnant scenario is not tied directly to the 
MO condition that cooling occur slightly before overlap. That 
condition is specific to the choice of cloud population and 
supernova rate used. (IX) argue that the cloud population will 
adjust to bring this about, but let us suppose that it does 
not.) More generally, then, the X-ray consistent conditions are 


10 


-3 


cm -3 ( 100 3/2 
TTT v P j 


J 51 



n = 


(28a) 



22 


R 


s 


228pc 


/ E 51 P 1/2 1 

MOOB j T 3/4 


(28b) 


1.23 [• 


E 51 P 
Q a 2 (100B) 3 


2/11 


(28c) 


In Figure 2, locus A is provided by equation (28a) with Tg » 1, or 
equation (13), and 100B “1. It provides the n Q (P) limit between 
allowed and overly bright remnant populations. 

The formal value of the areal filling factor at the cooling epoch 
is 


£ 1/2 p 13/3 

Qa c - 1-43 x 10 4 — rpTTo • (28d) 

(100B) 10/3 


Low and intermediate dilutions accompany c > 0.5, high dilution 
Qa, c ^ 0.5. If the M101 remnants follow an evaporative evolution 
yet provide B < 10~^, then unless P is small, remnants will 
overlap long before cooling. Very low densities and large radii 
are then the rule, in addition to remnant merger to form a hot 
corona. 

These conditions have been made already consistent with the 
X-ray brightness so the results seem to ascribe a maximum value of 
n Q consistent with P, for any value of P . If however, remnant 
overlap occurs while the temperature is still in the X-ray 
emitting regime (T > 5 x 10^K) , the corona is also subject to the 
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constraints of the McCammon and Sanders (1984) analysis. As we 
shall see, these evaporative evolutions leading to a hot corona 
are also too bright in X-rays, again unless the supernova power is 
less than normal. 

The conditions at the base of a galactic fountain or wind are 
like those of the source remnants when ~l/3 (Cox, 1985b). Thus 
the fountain temperature can be expected from equation (28c) to be 
approximately 

9 E S1 p5 2/11 

X 6 f - 1.23 [ — - 3 ] 2 / 1 1 (29) 

(100B) 

The McS surface brightness limits were given for 2o confidence so 

even though we are now considering halo emission we adopt their 
20 —2 

= 2 x 10 cm results to remain conservative. 


For < 10 ^K their surface brightness limit can be approximated 
by: 


-(8.55 

S < 10 


1*5/Tg f) _2 

’ erg cm z 


1 sr 1 


(30a) 


which corresponds to 


-(3.66-1.5/T 

P < 10 


6 , f 


) 


(30b) 



24 


The behavior of the exponential is sufficiently strong that unless 
E 5 J or 100B is very different from 1 consistency with equation 
(29) requires Tg f < 0.50, and 

P E 51 1/ ' 5 < 0.24 (IOOB) 3 / 5 . (31) 

The more precise criterion is shown in Figure 2 as the boundary 
between the allowed fountain and overly bright coronal regimes. 

The presence of B in equation (31) may seem extraneous since we 
have imposed the McS analysis. The X-ray result has in fact been 
used twice. The evolving remnants must not generate too many X- 
rays before overlap (introducing B) nor must the resulting 
fountain after overlap. It is represented in Figure 2 by 
intersection C. 

Consulting the now complete Figure 2 we find the highest 
allowed supernova power to be at intersection C in the 
intermediate dilution regime. Its numerical value is given by 
equation (31). Similarly, the highest power allowed to an MO 
state (critical dilution) is at intersection B, the numerical 
value provided by equation (26). High dilution possibilities have 
essentially this same upper limit to their power. Remnant 
density n Q — 0.1 can be accommodated by a further reduction in 
power by a factor of 2 to P ~ 0.05. Finally, low dilutions, 
leading to coronal temperatures in excess of 2 x 10 ^K and hot 
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winds with little X-ray emission, are still conceivable, but they 
lie off the upper left corner of Figure 2. Their exact 
delineation is outside the scope of this paper. 

Allowing a supernova rate which just barely prevents an X-ray 
luminous corona (as at intersection C in Figure 2) is a dangerous 
game. Owing to the temperature sensitivity of the emission, very 
small irregularities in local supernova rates, or cloud popula- 
tions for evaporation, could lead to local hot spots in the foun- 
tain that would violate the constraint. The MO state is thus a 
better bet as an upper limit on the rate. 

For this reason, and because it offers a specific standard 
for comparison, we finally present the parameters of the MO state 
which is just barely consistent with the present limit. It falls 
to future studies of the M101 remnant population to discover 
whether these results are an accurate portrayal. The parameters 
are: 


T gN R g 2 “ 330 yrs <15kpc) 2 E 51 29/26 (32a) 
R s - 70pc E 51 23/52 /T 6 3/4 (32b) 
n = 0.035 cm -3 T 6 5/4 /E 51 17/52 (32c) 
t = 1.52 x 10 5 yrs E 51 23/52 /T 6 5/4 


(32d) 
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£ * 1.34 x 10 37 erg s~* Egj 33732 /!^^ 4 (32e) 

f = 0.053 E 51 3/26 /T 6 3/2 (32f) 

Q a *= 0.0046 E 51 11 / 52 /T 6 11 / 4 = 1.01 f 11 / 6 . (32g) 

In this list, R s , n, £, and f are the shock radius, mean 
remnant density (hot phase only), age, luminosity, and radiated 
fraction of the total energy, as functions of the post shock 
temperature. Q A is the areal filling factor of the disk for 
remnants with post shock temperatures T and higher. Notice that 
the required supernova rate, at least that in the disk but outside 
superbubbles, is quite low, as are the remnant densities. Condi- 
tions at the cooling and merger epoch can be found from setting Q A 
= 0.5. 

V. OVERVIEW 

Given the propriety of the general model of Section II, the 
remaining assumptions acted to maximize the upper limit on the 
brightness parameter B. For the standard Sedov solution (b = 0) , 
the resulting limit on density within remnants is given by 
equation (15b). At the expected supernova rate, it marginally 
threatens the possibility that remnants evolve in an intercloud 
component of density 0.1 to 0.3 cm -3 . 
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For evaporative evolution of remnants, no consistent system 
exists at the expected supernova rate. The possibilities at lower 
power are summarized in Figure 2. Below normalized supernova 
power per unit area (Equation 12) P ~ 0.2, , a broad range of 
possibilities opens up. The critical dilution MO states and high 
dilution disk confinement of remnants are available for P < 0.09. 

Any further conclusions depend sensitively on one's 
disposition. If one favors MO states, one might be lead to 
believe E 5 J <,0.5, particularly after remnants have suffered cloud 
compression and cosmic ray acceleration losses. One is then 
driven to accept a low rate ( x sn > 150 years per disk area of 

radius 15kpc) for relevant supernova or to fault finding with 

the assumptions of Section II. Conversely one might be tempted to 
see yet another nail in the coffin of MO ISM models. Pushing then 
for return to a higher density intercloud component, one finds 
that the then maximum acceptable SN rate for b » 0 is only a 
factor of three above that which admits the MO state for b » 

—5/3. Perhaps that is comfortable but one must keep a wary eye on 
the conservative assumptions leading to the limit as well as 
future observations and refined analysis that will almost 
certainly push it downward. Those familiar with the disquieting 
X-ray properties of the LMC remnant population (e.g. Long, 1983; 
Cox 1985a) may just be further convinced that the remnants do not 
evolve as expected and that this type of study must wait until 
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individual remnants in the population can be studied. Finally, 
one might see further confirmation of the need to consider 
correlated SN explosions in clusters as an important component, 
possibly driving a hot non-radiative wind, reducing the residual 
SN rate with which the general disk must cope. 

While some useful constraints are placed by the Einstein 
observations, other potential conclusions are statistically 
marginal or just out of reach, and data of even this quality do 
not exist for any suitable galaxy except M101. The observational 
situation should improve tremendously in the foreseeable future, 
however. The Advanced X-ray Astronomy Facility (AXAF) will of 
course have greatly improved angular resolution and sensitivity, 
but on a shorter time scale Germany's Rontgen Satellite (ROSAT) 
should be very well matched to this particular task. Its 2 keV 
upper energy cutoff is no disadvantage, and optimization for lower 
energies provides a net effective area almost five times that of 
Einstein at 280 eV (and almost half that of AXAF). A much more 
important consideration is that at 20" half-power radius, the IPC 
angular resolution is a factor of three better than Einstein . 

This allows far more effective removal of point source 
contributions and permits observations of more distant galaxies, 
thus greatly increasing the sample size. 

The faster optics and improved IPC should result in an 
extraneous detector background which is usually quite negligible 
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compared to the diffuse X-rays. This and the more stable and 
linear operation of the new IPC will greatly reduce the systematic 
uncertainties which complicate analysis of diffuse data from 
Einstein . Another major benefit will be improved energy 
resolution which permits unambiguous separation of the .100-. 280 
keV X-rays from those in the .500-1.00 keV band. This will allow 
some useful temperature limits to be placed on the emitting 
material. 
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Table 1 


Einstein IPC Observations of M101 a 


Concentric Ring Number of Width 5 Arcmin 


Channel 


1 


2 


3 


4 


5 


6 

2 

11.9 

(3.5) 

10.2 

(1.1) 

10.5 

(0.7) 

12.6 

(1.0) 

10.7 

(1.7) 

00 

• 

Ln 

(1.1) 

3 

11.8 

(3.2) 

9.8 

(0.9) 

9.4 

(0.6) 

9.1 

(0.7) 

8.6 

(1.5) 

8.8 

(0.9) 

4 

11.3 

(2.7) 

6.8 

(0.7) 

6.6 

(0.5) 

6.7 

(0.6) 

7.1 

(1.2) 

5.7 

(0.7) 

5 

8.1 

(2.5) 

5.3 

(0.7) 

3.6 

(0.4) 

5.1 

(0.5) 

6.3 

(1.1) 

4.9 

(0.6) 

6 

3.9 

(2.1) 

4.2 

(0.6) 

3.2 

(0.4) 

3.2 

(0.5) 

3.0 

(0.9) 

3.0 

(0.6) 

7 

5.9 

(2.2) 

2.3 

(0.5) 

2.0 

(0.4) 

2.0 

(0.4) 

1.4 

(0.8) 

2.6 

(0.6) 


a (Rates In units of 10 -5 counts (arcmin^s) “1 are followed in brackets by 


their la uncertainties.) 
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FIGURE CAPTIONS 

Fig. 1 . Upper Limits to the Brightness Factor B. The limits are 
shown versus the density evolution factor b for various 
values of the hydrogen column density N H (in units of 
10 ^® cm -2 ) assumed to overlie the emitting material in 
M101. The upper four curves are for ring 1 (minus 4+5 
as background), the lower four for rings 2+3. 

Fig. 2 . Constraints on Populations of Evaporative Remnants in 

M101 for B < 0.01. P is the normalized supernova power 
per unit area of equation ( 12 ) while n Q is the average 
density of the hot component within a remnant when the 
temperature is 10^K. Solid lines are for = 0.5, 
while dashed are for E 5 J *= 1. Locus A is from equation 
(13). Intersection B is from equation (26). (The 
remmant properties for conditions at this intersection 
are provided by equation (32).) Intersection C is from 
equation (31). The boundary to acceptable fountains 
follows from combining equations (13) (29) and (30b). 

The MO state locus, from equation (23), separates the 
intermediate dilution regime of fountains and radiative 
coronae from the high dilution regime of disk-confined 


remnants . 




b 


Figure 1 
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